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ABSTRACT 

Aims. We study the relations between stellar kinematics and chemical abundances of a large sample of RAVE giants in search for 
selection criteria needed for disentangling different Galactic stellar populations, such as thin disc, thick disc and halo. A direct com- 
parison between the chemo-kinematic relations obtained with our medium spectroscopic resolution data and those obtained from a 
high-resolution sample is carried out with the aim of testing the robustness of the RAVE data. 

Methods. We select a sample of 2167 giant stars with signal-to-noise per spectral measurements above 75 from the RAVE chem- 
ical catalogue and follow the analysis performed by Gratton and colleagues on 150 subdwarf stars spectroscopically observed at 
high-resolution. We then use a larger sample of 9131 giants (with signal-to-noise above 60) to investigate the chemo-kinematical 
characteristics of our stars by grouping them into nine subsamples with common eccentricity (e) and maximum distance achieved 
above the Galactic plane (Z^.^,^). 

Results. The RAVE kinematical and chemical data proved to be reliable by reproducing the results by Gratton et al. obtained with 
high-resolution spectroscopic data. We successfully identified three stellar populations which could be associated with the Galactic 
thin disc, a dissipative component composed mostly of thick-disc stars, as well as a component comprised of halo stars (presence of 
debris stars cannot be excluded). Our analysis, based on the e-Z„,^x plane combined with additional orbital parameters and chemical 
information, provides an alternative way of identifying different populations of stars. In addition to extracting canonical thick- and 
thin-disc samples, we find a group of stars in the Galactic plane ( Z^ax < 1 kpc and 0.4 < e <0.6), which show homogeneous kine- 
matics but differ in their chemical properties. We interpret this as a clear sign that some of these stars have experienced the effects of 
heating and/or radial migration, which have modified their original orbits. The accretion origin of such stars cannot be excluded. 

Key words. Galaxy: abundances - Galaxy: evolution - Galaxy: structure - Galaxy: kinematics and dynamics 



1. Introduction 

The chemical enrichment of the Universe is one of the main 
thrusts of modern astrophysics, and the Milky Way (MW) can be 
seen as the Rosetta stone of this evolution. In particular. Galactic 
archeology, i.e., the combined study of kinematics and chemical 
composition of stars of different ages, has recently become one 
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of the cornerstones of research in galaxy formation. The main 
goal in this new field is to reconstruct the formation history of 
our Galaxy by analyzing its fossil chemical records and kine- 
matical information. An important sub-product of this kind of 
study is also to provide observational constraints to models of 
galaxy formation in general. While there has been considerable 
progress in the past years to reproduce realis tic dis c galaxies in 
cosmologial simulations (e.g. Guedes et al. 1201 iL Brook et al. 
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l2012aL Stinson et al. 1201 ih . the outcome heavily relies on the 
advanced feedback models involving a cons iderable number of 
free p arameters (Piontek & Steinmetz 1201 1[ Scannapieco et al. 
I20T 1|) which need to be calibrated using empirically determined 
relations. 

Two main paths have been taken in the last years in the 
field of Galactic archeology. On the one hand high-resolution, 
high signal-to-noise (S/N) spectra of small/medium stellar sam- 
ples (~100-1000 stars) have been used to unveil the chemo- 
dynamical properties of the different galactic com poneri t s (e.g., 
Gratton et al. [l996; 20QS l2003l Fuhr mann [1991 120081 
Adibeky an et al~ [201 U Kordopatis et al., 1201 Ibi Bensbv & 
Feltzing 120121 and references therein). Most of these studies 
have the disadvantage that they are based on pre-selected sam- 
ples defined according to strict kinematical criteria and hence 
suffering from biases that are difficult to quantify. 

On the other hand, large spectroscopic surveys, suc h as th e 
Geneva Copenhagen Survey (GCS, Nordstrom et al. |2004| - 
~ 16000 stars), the Sloan Extension for Galact ic Understanding 
and Exploratioifl (SEGUE, Yanny et al. l2009l) and the RAdial 
Ve locity Experiment (RAV E, Steinmetz at al. l2006l Zwitter et 
al. l2008L Siebert et al. l201 ih take advantage of their large sample 
statistics (hence sampling a large parameter space) to compen- 
sate for the much less precise measurements of the abundances 
and stellar parameters, which are typical for medium/low spec- 
tral resolution (or photometry, in the case of GCS). 

Most of the surveys quoted above aim at estimating the 
[a/Fe] ratios of large samples of stars, which, to a first ap- 
proximation, can be used as a proxy for the stellar relative age 
(this seems to be valid even when radial migration is present, 
see Sc honrich & Binney i2009bi Minchev, Chiappini & Martig 
l2012d) . The uncertainty in the [a/Fe] rati os com i ng from the 
SEGUE abundance pipeline (Lee et al., l2008at l2008bl) are 
around 0.2 dex, for S/N above ~ 20. For the GCS much less pre- 
cise estimates, based on Stromgr en ind ices as proxy for [a/Fe] 
were obtained (Casagrande et al.. l201 ll) . 

Although the original goal of the RAVE survey is to obtain 
radial velocities for up to a million stars, its spectral measure- 
ments turned out to be very useful for chemo-dynamic analy- 
sis. Recently, the first results of the RAVE abundance pipeline 
were published (the RAVE Chemical Catalogue - Boeche et al., 
1201 Ih . showing that it is possible to measure up to seven indi- 
vidual elements from spectra with S/N above 40, and at least 
three of them from spectra with S/N ~ 20. The RAVE-chemical 
pipeline uncertainties for [a/Fe] are similar to those of SEGUE 
- around 0.2 dex. 

A first study of the kinematic-abundance properties of the 
MW thin and thick discs using RAVE data has been carried 
out by Karatag & Klement (2012), using results coming from 
the main RAVE pipeline (not from the chemical pipeline). The 
authors used a sample of ~4000 main-sequence stars (with log 
g > 3) from the second RAVE data release and classified the 
Galactic disc populations according to their distribution on the 
Viot-[M/H] plane, using the so-called X stellar population param- 
eter defined by Schuster et al. (fl99 3fl Although the authors suc- 
cessfully reproduced the mean kinematic properties of the thin- 

' For recent results w ith SE GUE data see Cheng et al.. 120121 (-7000 
stars ), Schl esinger et al. l201 ll ( ~400() G- and 23000 K-dwarfs), Bovy 
et a l I2OI II (-30000), Lee et al. l2011l (-17000 stars) and Liu & van de 
Ven l2012l f~27500 stars). 

- This parameter reflects the fact that thin disc stars tend to be more 
metal-rich and faster rotating than stars in the thick disc, although it 
is clear that both components strongly overlap both in metallicity and 
Galactic rotation velocity. 



and thick-di sc pop ulations as compared to previous work (e.g., 
Veltz et al., l2008h . their results should be taken with caution. 
Indeed, the metallicities and in particular the [a/Fe] estimates 
coming from the RAVE main pipeline are not free of consider- 
abl e error s and systematical effects (see a discussion in Siebert et 
al.. l20Tlh . which might be a problem for a classification method 
relying on metallicity, as the one adopted in Karatag & Klement 
(S2012) . 

Recently, the RAV E project adopted a new pipeline 
(Kordopatis et al. 1201 ll and Kordopatis and the RAVE collab- 
orators, in preparation) to estimate the stellar parameters val- 
ues as effective temperature, gravity and metallicity (Tes, \ogg, 
[M/H]), which is free of most of the systematic errors cited be- 
fore. Such values are required as input to the RAVE chemical 
pipeline which estimates the chemical elemental abundances. 
To date, the RAVE-chemical catalogue constitutes the largest 
sample for which individual chemical abundances are available. 
Here, we exploit the RAVE chemical catalogue in the frame- 
work of abundance-kinematic correlations of stars in the Galaxy. 
We aim at disentangling the different Galactic stellar compo- 
nents, as well as at identifying stars which were accreted or 
experienced heating and/or radial migration along the evolu- 
tion of our Galaxy. Possible processes that have been studied 
by means of N-body simulations include stel lar dif fusion driven 
by transient spiral arms (Sellwood & Binnev l2002[ Roskar et al. 
2008), by the interaction between spiral arms and the Galactic 
bar (Minchev & Famaey 1201 Ol Minchev et al. 1201 lal Brunetti 
et al. I2OI ih or d epositi ng stars into the Gal actic d iscs by merg- 
ers (Abadi et al. l2003l Villalobos & HelmiHool Di Matteo et 
al. l20T2l) . Such a disc stining can give rise to a number of phe- 
nomena in the chemo-kinematical properties of a galaxy, such 
as flatt ening in radial metallici ty grad ients (e.g., Schon rich & 
Binney l2009al Minchev et al. 1201 lal Pilkington et al. l2012t) . 
exten ded stellar den sity pro files (e.g., Sanchez-Blazquez et al. 
l2009l Minchev et al. l2012al) . and can have profound impact on 
the way we interpret the abundance-kinematic coiTelations in 
our own Galaxy as recently discussed by Minchev, Chiappini 
& Martig (I2012d) . 

From the observational side, the main difficulty has been the 
proper disentanglement of the local thin and thick discs. No se- 
lection criterion is free from biases, be it kinematical or chemi- 
cal (some drawbacks of a separation based purely on the [a/Fe] 
ratios in the framework of the SEGUE sample are discussed in 
Brauer et al. 2012, in preparation). Recently Bovy et al (20131) 
even argued that there is no clear dichotomy between the thin 
and the thick disc, a resu h in c ontrast to previous findings in 
similar studies (Veltz et al. l2008l) . It thus seems that the best way 
to address this intricate problem is to study the properties of a 
sample spanning the largest possible parameter space, both in 
kinematics and in chemistry. 

In the present work we make such attempt by using the best 
quality data of the RAVE abundance catalogue. Our two primary 
goals are: a) to show that the RAVE data give results consistent 
with those found from higher-resolution, higher-S/N analysis for 
the chemo-dynamic al properties of nearby stars; and b) to prop- 
erly investigate the abundance-kinematical properties of the thin 
and thick discs with a sample 10 times larger, and which covers 
a much larger volume (from ~ 100 pc to ~ 3-4 k pc), th an the one 
adopted in the pioneering work of Gratton et al. (|2003|) . hereafter 
G03. 

The layout of the paper is as follows: in Section 2 we discuss 
our sample selection. Section 3 describes how the orbital param- 
eters were computed. In Section 4 we carry out the analysis of 
the kinematic-abundance properties of the thin and thick discs 
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divided according to the kinematic selection criteria adopted in 
G03 and show that the RAVE abundance pipeUne gives trust- 
worthy results, despite the medium resolution of our spectra. 
In Section 5 we propose a different approach to the problem of 
disentangling different Galactic populations. Finally, discussion 
and conclusions are presented in Section 6. 



2. Data and Sample Selection 

This work makes use of the full suite of the RAVE data products, 
namely radial velocities, stellar parameters, distance estimates 
and chemical abundances. Information on the radial velocities 
and proper motions come from Siebert et al. (1201 ll) . distances 
have been computed following the recipe described in Burnett 
& Binney ( 1201 ll) . 

Recently RAVE adopted a new pipeline (Kordopatis et al., 

1201 ih for estimating the stellar parameters values. This pi peline 
makes use of the codes M ATISS E (Recio-Blanco et al., |2006|) 
and DEGAS (Bijaoui et al.. l2012l) optimized for the Can triplet 
region. These stellar parameters are adopted for the present 
work and a more comprehensive description of the data will be 
presented in the next RAVE data release (Kordopatis and RAVE 
collaborators, in preparation). 

The present chemical abundances internal data release em- 
ployed for this study contains elemental abundances for 245649 
MW stars (for a detailed des criptio n of how these abundances 
are obtained, see Boeche et al. l2011l) . Roughly 10% of the RAVE 
stars have been observed more than once. In order to avoid us- 
ing more than one estimation per star, for the re-observed stars 
we here adopt only the values derived from the spectrum with 
highest S/N. 

For this first paper using the RAVE abundance catalogue, we 
select stars with the highest quality spectra and abundances. We 
additionally require that such stars cover a chemical-abundance 
range as large as possible (both in metallicity and [a/Fe] ra- 
tios). To satisfy the first condition, we select spectra on which 
the code MATISSE converged to a single point of the param- 
eter space {Tes,logg, [M/H]), having high signal-to-noise ratio 
(at least S/N=60), which are well fit by the reconstructed spec- 
tra of the chemical pipeline (x^ < 1000) and with no contin- 
uum defect^ (frac > 0.99). We also take care that the selected 
spectra are not peculiar (classified as normal stars accor ding t o 
the three classification flags described in Matijevic et al. l20I2h . 
To satisfy the second condition, we select only cool giants and 
avoid dwarf stars. Indeed, owing to their weaker lines, the abun- 
dance measurements in dwarfs are more uncertain, especially for 
[X/H] < -1.0 dex. Conversely, for cool giants, thanks to their 
intense absorption lines, the RAVE abundance pipeline is able 
to measure chemical elements down to [X/H] - -2.0 dex. We 
exclude giants with \ogg < 0.5, to avoid any possible effects 
due to the boundaries of the learning grid used for the auto- 
mated parametrization (see Kordopatis and RAVE collaborators, 
in preparation). 

In the present work we consider the chemical abundances of 
Fe and a elements (the latter obtained as average of [Mg/H] and 
[Si/H]). These elements have the most precise abundances of all 
the elements included in the catalogue and, thus, are the most 
suitable ones for comparisons with the G03 work. 



' The frac parameter described in Boeche et al. (1201 ih gives the frac- 
tion of pixels that are non-defective and represents the goodness of the 
continuum fitting. 
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Fig. 1. logg versus T^g values for the SN75 RAVE sample (see 
text) with overplotted isochrones of 3 Gyr (left panel) and 1 Gyr 
(right panel) with metallicities [M/H] =0.0 dex (solid lines) and 
[M/H] =-1.0 dex (dashed lines). Black points, grey crosses and 
open squares represent the thin disc, dissipative and accretion 
components as defined in Sec.|4] 



We decided to work with two samples of difl'erent S/N: the 
main one contains only spectra with S/N>75 (the "SN75 sam- 
ple") and the second one with S/N>60 (the "SN60 sample"), 
used when a better statistic is needed because it contains a larger 
number of starfl Both samples satisfy the following criteria: i) 
giant stars with gravity 0.5 < log^ < 3.5; ii) effective tempera- 
ture 4000 < reif(K)< 5500; iii) high quality data (spectra with 
X' < 1000 and frac > 0.99). With these criteria, the SN75 and 
the SN60 samples count 2167 and 9131 stars, respectively. 

These stellar parameter constraints are optimal for selecting 
cold giant stars from the RAVE sample and avoiding the hot gi- 
ant stars of the horizontal branch. Figure [Tlshows where our s am- 
ple lies on the reff-lo g g plane. Also plotted are the isochrones of 
Marigo et al. (I2008h for metallicities 0.0 and -1.0 dex and ages 
between 3 and 10 Gyr. Figure |2] shows the spatial distribution 
of our sample on the xcai-ycai (Galactic disc) and xcai-zcai (ver- 
tical direction) planes, where the Sun's location is indicated by 
the intersection of the dashed lines in each panel. 

3. Orbital parameters 

To be self-consistent with the new stellar parameters adopted 
in this work (from the new RAVE pipeline - Kordopatis et al. in 
prep.), we recomputed the chemic al abun dances using the RAVE 
chemical pipeline (Boeche et al. |20I lb and new distances for 
our st ars by using the method described in Burnett & Binney 
(|20I ih improved to take in account interstellar extinction and 
reddening (Binney et al., in preparation). We next computed 
the 3D space velocities u, v and w along the Galactic coordi- 
nates XGai, jGai and ZGai- We corrected the velocities for the lo- 
cal sta ndard of rest velocity as derived by Dehnen & Binney 
(Il998h . In order to obtain additional orbital parameters, we nu- 
merically integrated the orbits of stars by using the code NEMO 
(Teuben,[l99l), given their stellar distances and velocities. For 
the Gal actic p otential we adopted the model n.2 by Dehnen & 
Binney (Il998h . which assumes Ro=8.0 kpc, circular velocity at 
the solar circle Vf.(Ro)=217 km s and disc surface density 
S=52. 1 Mg pc"- (even when potentials n. 1, 3 or 4 by Dehnen & 
Binney ( 1998) were employed, the results did not change signif- 
icantly). From the integrated Galactic orbits we extracted useful 



We verified that all the results of this work are valid for both sam- 
ples. 



4 



C. Boeche et al.: The relation between chemical abundances and kinematics of the Galactic disc with RAVE 



"1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 
I 
I 




J I I I I I I I I I I I I I I I I I I I I I I I I L 



~i — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



□ 



J I I I I I I I I I I I I I I I I I I I I I I I I L. 



-10 



-8 

XGal(kpc) 



-6 



-10 



-8 

XGal(kpc) 



Fig. 2. Distribution of the SN75 RAVE sample in the xcai - ycd\ and xcai - Zcai planes (left and right panels, respectively). The 
dashed lines cross at the Sun's position. Symbols as in Figure[T] 



quantities, such as apocentre, pericentre, Rp, eccentricity, e. 
Galactic rotation velocitjO Viot, and the maximum vertical am- 
plitude, Zn,ax- Ra and Rp are the maximum and minimum dis- 
tances from the Galactic centre a given star obtains during a 
revolution of In radians around the Galactic centre, measured 
on the Galactic plane. Similarly, Zmax is the maximum altitude 
reached by a star along its orbit. The eccentricity is defined as 
e^(R„-Rp)/{R„+Rp). 

In order to estimate errors in the orbital parameters we 
followed the G03's method, i.e., we performed a Monte Carlo 
simulation: for each star we computed 100 times the orbital 
parameters by changing every time the Galactic coordinates 
■^Gab ycah ZGai and the velocities u, v, w, according to their esti- 
mated errors. From these 100 orbits we computed the standard 
deviation for each orbital parameter. Despite the large distance 
for some of our stars (up to 3 kpc, see Figure |2] and Figure |5] 
top-right panel) for which the proper motions (and therefore the 
tangential velocities) are very uncertain, the orbital parameters 
show reasonably small variations when errors are taken in 
account. As showed in Figure [3] the errors in eccentricity are 
smaller than 0.2 for most stars, while for Rp they are smaller 
than 1 kpc. Similar behaviour is seen in the other parameters. 
With such errors we are confident that we can use these orbital 
parameters to discriminate between thin- and thick-disc stars, 
with only moderate contamination. 

During the preparation of this manuscript we used four dif- 
ferent versions of distance estimations (two prelimin ary an d to 
definitive, following the methods of Zwitter et al. 1201 Ol and 
Binney et al. in preparation) an d two d ifferent chemical catalog 
versions (DR3 by Siebert et al. 1201 ll and DR4 by Kordopatis 
et al. in preparation) for our analysis. The results of the present 




Fig. 3. Estimated uncertainties versus their respective stellar or- 
bital parameters value for the SN75 sample. Dashed, dotted and 
dash-dotted lines represent 100%, 50% and 25% errors, respec- 
tively. 



^ Following the common use, the name "Galactic rotation velocity" 
refers here to the azimuthal velocity in a cilindrical coordinate system, 
computed as V|ot = • ^ . 



work (see next sections) hold for all the data set we employed, 
confirming the robustness of the selected sample and the results 
obtained. 
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4. The kinematical criteria of G03 applied to the 
RAVE data 

In G03 the authors used 150 field subdwarfs and early subgiants 
with accurate parallaxes and kinematics and divided them into 
three subsamples according to pure kinematic criteria. With the 
purpose of validating the RAVE abundance pipeline, as well as 
our orbital parameters, we follow here the analysis done in G03 
and divide the RAVE samples into three subsamples according 
to the following criteria; 

1. The thin disc component including stars whose orbits 
have low eccentricity and low maximum altitude from the 
Galactic plane, i.e. e < 0.25 and Z^ax < 0.8 kpc. It will be 
represented as grey crosses in the figures that follow. 

2. The dissipative collapse component consisting of stars with 
Vrot > 40 km s"' . Stars belonging to the thin disc component 
just defined are excluded. This sample includes part of the 
thick disc, as well as of the halo. Through all the figures this 
sub-sample will be represented with black dots; 

3. The accretion component of non-rotating or counter-rotating 
stars. These stars satisfy the constraint Vi-ot < 40 km s"'. 
This population is very likely composed both from halo stars 
and accreted debris (i.e. stars that do not share the rotation of 
the disc components). It is represented in the figures as open 
squares. 

The criteria mentioned above deviate somewhat from the 
criteria used in G03, who defined the thin disc component via 

V-^max + 4e2 < 0.35 and had an additional constraint on the 
apogalactic radius Ra < 15 kpc for the dissipative compo- 
nent. The thin disc stars selected by this constraint cannot have 
e > 0.175 if they lie on the Galactic plane and they must have cir- 
cular orbits if their Z,nax = 0.35 kpc. Indeed, while the above cri- 
teria led to differences in the chemical composition that likely re- 
flect real differences in the stellar populations (as will be shown 
in the next sections), the detailed specification is to some extent 
arbitrary, as already pointed out by G03. Our modification al- 
lows us to generalize the selection criteria in the e - Zmax plane, 
as we will show in the next section. Furthermore, our sample 
covers a considerably larger volume than the Hipparcos sphere 
probed by G03, resulting in a relatively low number of stars with 
near-circular and coplanar orbits when the G03 original selection 
criteria is employed. 

Indeed, the G03 sample is composed of dwarf stars, cov- 
ering a small spatial volume (< 100 pc from the Sun) in or- 
der to include only stars with accurate parallaxes. Moreover, 
their sample was drawn from the Hipparcos catalogue for which 
metal-poor stars were preferentially selected (with [Fe/H] <- 
0.8 dex). Hence, this sample includes a rather large number of 
high-proper-motion stars, resulting in a strong kinematical bias 
favouring objects on highly eccentric orbits, with low Galactic 
rotation velocities, and with either large apogalactic or small 
perigalactic distances. 

The RAVE sample considered here, on the other hand, is 
composed by giant stars, thus, probing a large volume of space 
(up to 3 kpc from the Sun). Most importantly, the RAVE sample 
is rich in high metallicity (> -1.0 dex) stars. Moreover, while 
the G03 sample is composed of different high-resolution spec- 
troscopic subsamples, partly obtained by the authors and partly 
found in the literature, the RAVE sample adopted here is very ho- 
mogeneous. Finally, our samples SN75 and SN60 have respec- 
tively ~14 and ~60 times the size of the G03 sample. 

In Figure |4] we show number density contours of the max- 
imum height achieved above the Galactic plane, Z^ax versus 
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Fig. 4. Density distribution on the e - Zmax plane for our sam- 
ple SN75 sample, with 2152 stars. The dash-dotted white curve 
shows the G03 thin-disc selection criteria, whereas the dashed- 
dotted black line frames the area defined by our modification to 
the e and Z^ax constraints. The contours contain 34%, 68% and 
95% of the sample. 



the orbital eccentricities, e for the sample SN75. The criterion 
V-^max + 4e- < 0.35 adopted by G03 to define the thin-disc sam- 
ple is represented by the quarter of an ellipse on the e-Z^ax plane 
(white dash-dotted curve), the thin disc defined by our criterion 
is shown as black dash-dotted line. 

With the above prescriptions, the thin-disc component con- 
sists of 1079 stars, the dissipative component of 1024 stars, and 
the accretion component of 64 stars. The results of this division 
are shown in Figures [1] |2] E] |6] |7] and |8] The thin disc, dissi- 
pative and accretion components show different chemical signa- 
tures, similar to the ones seen in G03, with the difference that 
our sample extends up to solar abundances. 

Figure |5] shows the main properties of the different samples. 
While the accretion component appear scattered over a wide area 
of the Toomre diagram (and it would extend to higher v, since 
such non-rotating population must hold rotating and counter- 
rotating stars in equal amount) the thin disc and dissipative com- 
ponents clump and overlap on a smaller area. With respect to the 
distances, the outcome of our sample selection is that the thin 
disc stars cover a smaller volume than the dissipative sample, 
whereas the accreted stars prefer larger distances (between and 
3 kpc) from the Sun. The eccentricity distribution of each sam- 
ple is also shown. It is interesting to notice that the eccentricity 
and distance distributions found here are qualitatively similar to 
those found with the SEGUE G-dwarf sample used in Lee et al. 
(2011), even though the latter authors adopted a pure chemical 
criterion to divide their sample into thin- and thick-disc stars and 
even though our distance and eccentricity distribution is likely to 
be affected by our selection criteria. 

Indeed, the RAVE giant sample covers essentially the same 
volume as the SEGUE dwarf sample (SteinmetZ i2012i) . The sim- 
ilarities between our thin and dissipative sa mples and the thin 
and thick disc ones obtained by Lee et al. (1201 ih . is by itself 
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Fig. 6. Perigalacticon, Rp, eccentricity, e, apogalacticon, Ra, and 
Galactic rotation velocity, Vrot as a function of [Fe/H] for the 
three subsamples. Symbols are as in Figure |5] This figure corre- 
sponds to Figure 4 of G03. 



reassuring. Indeed, two completely different surveys, using dif- 
ferent tracers (giants vs. dwarfs), with chemistry, distances and 
orbital parameters computed by different pipelines, analyzed in 
completely different ways (kinematically selected vs. chemically 
selected), still give very similar answers for the mix of Galactic 
stellar populations within ~2-3kpc from the Sun. 

In Figure |5] we also show the [cc/Fe] distributions for each 
of our samples. An offset of ~0. 1 dex is seen between the thin 
and dissipative samples, with however a considerable overlap 
(as expected even according to pure chemical evolution mod- 
els, where the difference in [a/Fe] between the t hick an d thin 
discs is a function of metallicity - see Chiappini |2009|) . Note 
that the absolute values of the mean [a/Fe] ratios of the thin 
disc and thick disc (dissipative) distributions are in good agree- 
ment with the mean values reported by high-resolution studies 
(Bensby and Feltzing [20121 and references therein). These val- 
ues ar e, how ever, lower than the mean values obtained by Lee 
et al. (1201 lb with the SEGUE sample, namely: [a/Fe] -0.1 for 
the thin disc, and ~0.35 for the thick disc (a detailed comparison 
between RAVE and SEGUE results is beyond the scope of the 
present paper, and will be presented by Brauer et al. in prepara- 
tion). 

In Figures |6] |7] and |8] we successfully reproduced G03's 
Figures 4, 5 and 6, respectively, with some additional features. 
The RAVE data support some of the G03 results, such as the 
existence of correlations between [Fe/H] and Rp, e, Vmi, and 
the apparent absence of a correlation with Rg (Figure |6]l. The 
dissipative component (black dots) has a wide range of eccen- 
tricity, and moderately high abundances. Also seen in Figures |6] 
and |7] is an apparent drop in the black point density around 

[Fe/H] 1.0 dex and [a/H] 0.7 dex (see also Figure|9]and 

Figure \W\ top and middle panels). The standard deviation of 
the black points is ctj^^h] =0.2 dex (after the rejection of the 
accretion component stars polluting the dissipative component 
at [a/H] <-1.0), which is very close to the error in abundance 



Fig. 7. As in Figure |6] This figure corresponds to Figure 5 of 
G03. 
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expected for [a/H] (0.17dex). This support the weak correlation 
shown with Viot and -0.7 dex may represent the lower limit of 
the thick disk ff-abundance. 

In Figures |9] and [To] we show the trend of [Fe/o'] with [a/H], 
and [a/Fe] with [Fe/H] (similar to Figures 2 and 3 of G03, re- 
spectively). 

Both the thin and dissipative components show a decreasing 
[a/Fe] abundance ratio with increasing metallicity. However, the 
thin-disc component shows systematically larger [Fe/a] ratios 
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Fig. 5. Toomre diagram (upper left) for each of our samples: thin disc stars (grey points), dissipative component (mostly thick disc 
stars - black dots) and accreted component (squares). The light grey semi-circle indicates the constant peculiar velocity (u^ + + 
vi;2)i/2 _ 4Q i^jjj g-i Also shown are the distance (upper right), [a/Fe] (bottom left) and eccentricity (bottom right) distributions of 
the thin disc (grey line) dissipative (solid black line) and accreted (black dashed line) components. 



than the dissipative one for [a/H]< -0.2 dex. This can be seen 
by the solid red line, which indicates the average of the black 
points in bins with variable size in [Fe/H] so that every bin con- 
tains 50 points. By plotting this fiducial in each panel, the differ- 
ence between the thin-disc and dissipative component is visible 
in Figure|9] For [a'/H]> -0.2 dex the thin and dissipative compo- 
nents seem to be chemically indistinguishable. Note that for the 
accreted component the abundance ratio remains essentially flat, 
and more importantly, shows systematically lower [a/Fe] than 
the dissipative component (see bottom panel of Figure [TOli. This 
result is not only in agreement with the original findings of G03, 
but is also similar to the recent high-resolution results by Nissen 
& Schuster (i20ia) . 



It is clear that although the G03 criteria lead to clear sub- 
populations, overlaps between the samples still exist. The a- 
enhanced stars at [Fe/H] < -0.7 dex of the thin-disc component 
(Figure \W\ top panel) probably belong to the dissipative com- 
ponent. On the other hand, the stars at [Fe/H]< -1.0 dex of the 
dissipative component are very likely accretion component stars. 
In fact, because the accretion component stars are non-rotating 
(they can be identified as halo stars by looking at the locus they 
occupy in the Toomre diagram, top-left panel of |5) their average 
Vrot must be zero. By mirroring the square points with respect 
to y,.ot = (v = -217 km s"'), we can infer that many black 
points (dissipative component) could also be considered accre- 
tion component stars. 



Despite of the small overlaps described above, the three 
groups of stars assigned to the different components via kine- 
matical criteria, show different chemistry: 

- thin disc component: the distributions in Ra, Rp, e and Vrot 
span a small range of values because they are limited by the 
criteria e < 0.25 (nearly circular orbits) and Z^ax < 0.8 kpc 
(close to the Galactic plane), i.e., stars are contained in a lo- 
cal volume with R - ^ + 1 kpc, see Figure|2] This component 
has an [ff/Fe] peak slightly above zero and tend to be more 
shifted to higher metallicities (figures |6] and |7). 

- dissipative component: for this subsample Vrot is downward 
limited by the constraint Vjot > 40 km s"' and stars which 
happen to have circular orbits (e < 0.25) are missing from 
this component in favour to the thin disc component. This 
creates a small bias against the dissipative component which 
diminishes the number of stars at low eccentricities, but does 
not completely remove them, since the orbital parameters 
can still span a wide range of values for Rp and e as well. 
On the other hand, stars at [Fe/H]< -1.0 dex are likely to be 
members of the accretion component. The dissipative com- 
ponent has more stars at high [a/Fe] with respect to the 
thin-disc component. It also shows a correlation between 
the abundances [Fe/H], [a/H] and the orbital parameters, as 
found also by G03. 

- accretion component: the kinematical criteria adopted for 
this sample imply that such objects have mainly high ec- 
centricity (Figure |5). This component covers a wide range 



8 



C. Boeche et al.: The relation between chemical abundances and kinematics of the Galactic disc with RAVE 




Fig. 9. Abundance ratio [Fe/o-] versus the abundance [a/H] for 
the thin disc component (grey crosses, top panel), the dissipa- 
tive component (black points, middle panel), and the accretion 
component (open squares, bottom panel), selected by using our 
modified criteria (see text). The red line represents the average 
[Fe/o'] of the dissipative component obtained by averaging bins 
of 50 points each, and the eiTor bars represent their standard de- 
viation. The red line is reproduced in each panel as a fiducial 
line. This figure corresponds to Figure 2 of G03. 



Fig. 10. As in Figure|9]but for the abundance ratio [a/Fe] now as 
function of [Fe/H]. This figure correspond to Figure 3 of GOB. 



in [Fe/H], including objects with metallicities clearly larger 
than the upper limit halo stars metallicity of ~ - 1 . 

The successful reproduction of the main results of G03 val- 
idates the kinematical and chemical data of the RAVE survey 
(despite the fact that they were obtained from medium-resolution 
spectroscopy) and allows us to push further our analysis. 

However, any selection criteria aiming at disentangling the 
thin and thick discs will suffer from the fact that these two com- 
ponents overlap in almost all parameters. Pure kinematical se- 
lection criteria of thin and thick disc return samples whic h par- 
tially overlap in che mical abundances (Bensby et al. 20031 120051 
Reddy et al. l2006l among others), whereas pure chemical se- 
lection criteria return sample s which p artially overlap in kine- 
matics (Navarro et al. 1201 ll Lee et al. 1201 ll) . Even if a clear 
kinematical separation between these two components did ex- 
ist in the past, it could have been heavily blurred by a number 
of agents. In the next Section we try an alternative approach 
with the aim of avoiding strict selection criteria, which hopefully 
can be more successful in providing more robust constraints to 
chemo-dynamical models. 

5. The e - Zmax plane 

As discussed in Section 1, there exist two general processes, re- 
sulting from the disc secular evolution, which could partially de- 
stroy the kinematical borders between the thin and thick discs: 
(1) kinematical heating (increase of stellar velocity dispersion) 
with time and (2) radial migration (change in the angular mo- 
mentum of stars and, thus, in their guiding radii). Another pos- 
sibility is (3) the deposition of material by accretion, making the 
situation even more complex since mergers also can give rise to 
(1) and (2) above. Processes (1), (2) and (3) will lead to different 
signatures: 



(1) Heating by transient ( Carlber g & Sellwood. 1 19851) and 
multiple (Minchev & Quillen, '2006*) spiral density waves, the 
Galactic bar (Minc hev & Famaey|2010j), giant molecular clouds 
(Jenk ins & Binnev ll990l) . and minor merger events (Quinn et al. 
Il993h would most probably increas^ the velocity dispersions of 
all three components (U, V, and W). An increase in the vertical 
velocity dispersion would then result in overlapping an initially 
cooler stellar population (such as the thin disc), for which the 
heating is more effective due to the cold orbits, with an initially 
hotter stellar sample (such as an old thick disc), which would 
not be affected much by the heating agents. Note that the above 
would be true irr especti ve of whether the thick disc were born 
hot (Forbes et al. 1201 ih or were preheated by mergers at high 
redshift (e.g., Villalobos and Helmi |2008|) by virtue of the thin 
disc being younger than the thick disc and the expected decrease 
of merger activity with redshift. 

(2) In contrast, as a star migrates (gains or loses angular 
momentum), information about its birth radius is lost and its 
kinematics can be virtually indistinguishab le from those of stars 
bom at the ne w radius (Sellwood & Binnev l2002i Minchev et al. 
l20TTbll20I2bh . In recent years radial migration (or mixing) has 
been recognized as an important process affecting galactic discs. 
Several radial migration mechanisms have been described in the 
literature: the effect of the corotation resonance of transie nt spi- 
ral density waves (Sellwood (fe Binnev l2002l Roskar et al. 120081 
Schonrich & Binney l2009al) . the effect of the non-linea r cou- 
pling between multiple spiral wave s (Minchev & Quillen. l2006h 
or ba r and spirals (Minchev et al. 1201 Ol 1201 lal Br unetti et al. 
1201% and the effect of minor mergers (Quillen et al. |2009l Bird 
etaLMOb- 



* Note, however, that this effect still needs to be better quantified for 
a MW-like galaxy. 
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(3) Deposition of material by accretion has been shown to 
comm only occur in cosmological simulations (e.g., Abadi et al. 

In this scenario tidal debris of satellites with orbits coin- 
ciding with the plane of the host disc can populate a disc com- 
ponent and, thus, blur the borders between the preexisting thin 
and thick disc. 

Such diffusion mechanisms make ineffective any kinematic 
criteria aiming at separating the thin from the thick disc, even in 
the unrealistic case of no errors in the kinematical data avail- 
able. The chemical properties acquired at birth, on the other 
hand, must be preserved. However, here the difficulty is that 
due to the shallow abundance gradient in the galactic disc, as 
well as the small range in the [a/Fe] variation among differ- 
ent galactic populations (at most 0.5 dex), the chemical differ- 
ences are subtle (unless one uses abundance ratios involving 
other chemical elements which present larger variations). In ad- 
dition, even if the differences in the star formation histories of the 
thick and thin d iscs would have led to chemical differences (see 
Chiappini l2009l) . large overlaps would still exist. This would re- 
main true if the disk was composed of several mono-abundance 
sub-components, as suggested by Bovy et al. (201 1 ). For these 
reasons, we here abandon the "selection criteria" approach and 
try to have a different view of the problem. 

Our previous analysis considered the distributions of the 
orbital parameters (such as Z,nax, Rp, e and Vi-ot) separately. 
However, these distributions are slices in a bigger chemo- 
kinematical, multidimensional space, in which the stars lie. It is 
therefore possible that some information could be missed when 
only the shapes of these distributions are considered, resulting in 
a mixup in our separation of different Galactic components. For 
instance, consider two stars with Vmt ~22Q km s"' ; one can have 
a circular orbit on the Galactic plane at the Sun radius Rq, and 
the other can have Rp < Rq, an eccentric orbit and vertical ve- 
locity V- large enough to reach Z^ax > 2 kpc. It could be hardly 
thought that such stars belong to the same population. To avoid 
such a trap, we further analyse our sample by grouping the stars 
by similar orbits using the e - Zmax plane. The eccentricity gives 
the shape of the orbits, whereas Z^ax tells us about the oscillation 
of the star perpendicularly to the Galactic plane. 

In Figure [TT] we divide the e - Z^ax plane in nine groups of 
stars, as indicated by the dashed lines, and we label them from 
(a) to (i). We have neglected stars with e > 0.6 because here 
we focus on the Galactic disc. In this way we have sorted the 
stars into "classes" of orbits: moving rightwards the eccentricity 
grows, moving upwards Z^ax grows (hence the vertical v~ veloc- 
ity increases). 

By dividing the e - Z^ax plane into nine regions, we have 
obtained stellar subsamples with narrow ranges in orbital pa- 
rameters. For each group we now plot the [a/Fe] versus [Fe/H] 
relation and the distributions of [Fe/H], [c/Fe], Rp, R,„, and Vrot, 
in Figures [T2]and[T3] respectively. Panels (a), (c) and (g) of the 
aforementioned figures give particularly interesting insights, on 
which we now focus. 

5.1. Identification of tliin/tliick/diffused stars 

In the following analysis we use the SN60 sample, in order to 
have a better statisticO- We first study the distributions of sub- 
sample (a), as defined in Figure[TT| The stars in this group (with 
eccentricities < 0.2 and low vertical velocities) show the ex- 
pected properties of a sample dominated by local thin-disc stars, 

' We verified that the results found in this work are valid and consis- 
tent for both SN75 and SN60 samples. 
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Fig. 11. The e - Zmax plane divided in nine panels labelled from 
(a) to (i). Here we use the SN60 sample, which counts 9131 stars. 



namely: a Fe-distribution peak at ~ -0.25 dex, a Vi-ot peak at 
~220 kms"', and a mean galactocentric distance/?,,, of ~7.5 kpc. 
The apparent low abundance of the Fe-distribution (which peaks 
at ~ -0.25 dex instead o f ~ - 0.05 dex of the local thin disk 
found by Casagrande et al. l20I Ih is due to the spatial distribution 
of the RAVE sample, which lacks of nearby stars and favours 
stars lying at zcai > 300 pc, where the metallicity distribution 
funct ion is shifted to lower values (also found by Schlesinger 
et al. l20Ill and consistent with th e predic tions of the chemody- 
namical model of Minchev et al. Nonetheless, the Fe- 

distribution in panel (c) is more metal rich in comparison with 
the subsamples of the other panels. 

By moving up from panel (a) to panel (g) the mean Rp and 
R,„ do not change significantly but their distributions become 
broader (see Figure [T3]top panel), whereas the Vrot distribution 
slightly decreases its mean from 220 km s"' to 200 km s"' (see 
Figure [12] bottom panel). The Fe-distribution is found to shift 
to lower abundances (with a peak around -0.6 dex. Figure [12] 
top panel), and the average [a/Fe] increases by ~0. 1 dex. All the 
above properties are indicative of a sample dominated by local 
thick disc stars (panel g). However, the presence of kinematically 
heated old thin disc stars cannot be discarded. Stars with high 
Zmax and low eccentricities must be on orbits which strongly os- 
cillate through the Galactic plane, suggesting that they have ex- 
perience d som e perturbations during their lifetimes. Feltzing & 
Bensby (l2008h proposed the same interpretation for a subsam- 
ple of 32 stars which would lie in panel (d) of our Figure [TT] 
The cause for this might be identified in the kinematical heating 
mechanism proposed by some authors (e.g. Villalob os & Helmi 
.2008;. Bournaud et al. 2009, Minchev et al. l20I2d for different 
processes) in order to explain the thickness of the disc. 

Focusing now on panel (c) of Figures [T2l and [T3](upper pan- 
els), we identify a population with high eccentricities but con- 
fined close to the disc plane (Zmax < 1 kpc, 0.4 > e > 0.6). 
Stars in this group have small perigalactic radii (Rp ~3 kpc) 
and small mean radii {R,„ ~6 kpc). Due to their high eccentric- 
ities and the large fraction of stars coming from the inner disc. 
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Fig. 12. Upper panel: Relative abundance [a/Fe] versus [Fe/H] for the stellar samples defined 
by panels (a) through (i) in Figure [TT] The histograms represent the Fe distributions with 
relative scales. Lower panel; Distributions of abundance [a/Fe] for the stellar samples defined 
by panels (a) through (i) in Figure [TT] The distributions are normalized over the total number 
of points contained in each panel (Ntot). 
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Fig. 13. Upper panel: Perigalacticon (grey lines) and mean radius (black lines) distributions 
for the stellar samples defined by panels (a) through (i) in Figure[TT] The dotted line indicates 
^=8 kpc. Lower panel Rotational velocities Vrot for the stellar samples defined by panels 
(a)-(i) in Figure [TTI The dotted line indicates Vrot=220 km s"'. 
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the observed Viot distribution (Figure [13] lower panel) peaks at 
~120 km s"', and has a tail extending up to more than 200 km 
s"'. In addition, stars in this panel show a broad Fe-distribution 
(Figure [12] upper panel) with a hint of bimodality. The latter 
differs significantly from the metallicity distributions in panels 
(a) or (b) (a Kolmogorov-Smirnov test showed that the proba- 
bility that the distribution in panel (c) is drawn from the popu- 
lations of panel (a) or (b) is both lower than 10"'^). Note also 
that, by moving upward from panel (c) to panel (i), the high 
[Fe/H] tail of the distribution progressively disappear, leaving 
a Fe-distribution qualitatively similar to the thick disc one. This 
means that in panel (c) there are most probably more than one 
population, i.e. stars that might belong to both the thin and thick 
discs. While stars with low Fe abundance ([Fe/H]< -0.5 dex) 
have all the characteristic to be identified as thick disc, stars with 
high Fe abundance ([Fe/H]> -0.4 dex) are likely to be thin-disc 
stars scattered outward from the inner part of the Galaxy. Being 
such stars drawn from panel (c) their kinematic is determined 
and thick-disc-like, but their chemical abundance is typical of 
the thin-disc. Therefore, they have no clear thin- thick-disc clas- 
sification. This feature is highlighted in Figure [14] where the dis- 
tributions in Vi-ot, Rm and [a/Fe] of the two tails are shown sep- 
arately. Kinematically, the stars belonging to the two tails show 
no significant differences, whereas they have distinct chemical 
abundances in [Fe/H] and [(r/Fe]. 

Such stars having thick-disc kinematics and thin-disc chemi- 
cal abundances might have been kinematically heated and/or mi- 
grated by a mechanism which scatters out stars from the inner 
parts of the Galaxy, and could be identified as t he grav itational 
actio ns of t he spiral arms (Schonrich & Binney l2009a[ Roskar 
et al. 120081) a nd/or th e Galactic bar (Minc hev & Famaev l20Toi 
Minchev et al. l2011al Brunetti et al. l2011l) . 

The identification of two populations coming from the in- 
ner parts of the Galaxy, lying in the Galactic plane, but having 
different metallicities, suggests that thick-disc and dynamically 
heated and/or migrated stars are separate populations with dif- 
ferent origin and evolution. This obse rvatio n is consistent with 
the pre vious findings by Wilso n et al. (1201 ih . Liu & van de Ven 
(1201 2h and Kordopatis et al. (1201 Ibl) . based on the eccentric- 
ity distribution of thick disc stars and challenges scenarios in 
which the thick disc formed solely by the outward migration of 
stars born in the inner disc (Schonrich & Binney, 2009a), while 
on the other hand supports the predic ted ex istence and action 
of suc h a mechanism (see also Grenon [l999l and Trevisan et al. 
1201 Ih . Our conclus ion is a lso in agreement with the recent work 
by Minchev et al. (1201 2cl) . who suggested that both mergers at 
early phases and the effect of a central bar at later times are nec- 
essary to explain the presence of stars with thick-disc chemistry 
and kinematics currently found in the solar vicinity. 

The fact that we find such stars in panel (c) does not ex- 
clude the possibility that heated or migrated stars might be found 
in other regions of the e - Z^ax plane. For instance, as men- 
tioned in Sec.[5]migrated stars can move away from their origi- 
nal Galactic radius and conserve their eccentricity, making them 
indistinguishable from the locally bom stars. In such case only 
the chemical composition can reveal the difference between mi- 
grated and local stars, but in our sample, the identification of 
such stars is challenging, because to draw a clear chemical sig- 
nature requires a number of measured elements and a precision 
in abundance higher than the ones provided by RAVE. 
On the other hand, kinematically heated stars would have high 
probability to get into high eccentric orbits. In doing so, they 
would enrich the poorly populated tail of the eccentricity distri- 
bution of the local stars, and broad the [Fe/H] distribution, mak- 
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Fig. 14. Distributions in Vrot and R,„ of the high metallicity tail 
([Fe/H> -0.4 dex, black fine) and low metallicity tail ([Fe/H]< 
-0.5 dex, grey line) of the distribution in Figure [12] top, panel 
(c). 



ing more likely to see a bimodality within the error This make 
panel (c) a more favourable region where to look for such stars. 

6. Discussion and conclusions 

We selected two samples of 2167 and 9131 giant stars from the 
RAVE chemical catalogue having good quality spectra (S/N>75 
and S/N>60, respectively called "SN75" and "SN60" sample) 
and known kinematics. We app lied to the SN75 the kinematic 
criteria by Gratton et al. (|2003|) (G03) aiming at disentangling 
different populations of the Milky Way. Following G03, we di- 
vided the SN75 sample into thin disc, dissipative component and 
accretion component. We veiified first that our data (obtained 
from medium resolution spectra) matches well the results by 
G03 obtained from high-resolution spectra, confirming the re- 
liability of the RAVE kinematic and chemical abundances. 

Our thin disc and dissipative (mostly thick disc) samples se- 
lected by a pure kinematical criteria turned out to have similar 
dista nce an d eccentricity distiibutions as those reported by Lee 
et al. 1201 ll for a G-dwarf SEGUE sample where the thin and 
thick disc stars were selected based on a pure chemical criterion. 

We also notice that thin disc, dissipation and accretion com- 
ponents partially overlap in several parameters. Such overlaps 
are to be expected if the processes currently debated in the liter- 
ature (accretion, heating and stellar migration) are at play during 
the evolution of the Milky Way. This is particularly true for thin 
and thick disc (the latter is identified with the dissipation com- 
ponent) which kinematically overlap one another in a way that 
makes it difficult (if not impossible) to find selection criteria ca- 
pable of disentangling them. This realization pushed us to drop 
the searching for better selection criteria and led us to look for 
an alternative approach. 

The novelty of this work is the introduction of the analysis of 
the e - Zmax plane[3 which permits us to group stars having sim- 
ilar orbits. In fact, eccentricity determines the shape of the orbit 
and Zmax determines the oscillation amplitude of the stars in the 
Galactic plane, i.e., the probability of the star to populate the re- 
gion typically occupied by the thick disc. We applied this analy- 
sis to the SN60 sample, and divided it in 9 groups on the e -Zmax 
plane. By studying the distributions of R„,, Rp, Vmt, [Fe/H] and 
[ff/Fe] of these groups, we found stellar samples which identify 
nicely the thin and thick discs, as well as a third sample lacking 
a clear thin- thick-disc classification. 

In particular, we identified an interesting subsample of stars 
with large eccentricities (0.4 < e < 0.6), low Z^ax (below 1 kpc), 
with guiding radii in the inner disc {Rp ~ 3 kpc and R,„ ~ 6 kpc). 



A diffe rent u se of the e-Z,„ax plane was previously done by Feltzing 
& Bensby l l2008h . 
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and which shows broad distributions in [Fe/H] and [a/Fe]. The 
[Fe/H]-poor tail of the distribution is composed of stars which 
have properties similar to thick disc stars, whereas the high 
[Fe/H] tail shares with them the kinematical signature but differs 
clearly in Fe abundance as well as in a enhancement {[a/Fe]). 
A sample of stars with similar orbits and different metallicities 
suggests the existence of a heating/migration mechanism which 
pushes stars from inner part of the Galaxy outward. However, we 
cannot discard the possibility that such stars have been kinemat- 
ically heated by merging events or be long to a merged satellite 
themselves (see e.g. Helmi et al. l2006l) . 

Our results support a number of previous works, which have 
shown that the thin and thick discs overlap kinematically. Even 
if a clear separation did exist in t he pas t, heating from spiral 
arms (e.g., Carlberg & Sellwood 1 19851 Minche v & Q uillen 
I2OO6I) . a central bar (e.g., Minchev & Fa maey 1201 Ol) . giant 
molecular clouds (e.g., Jenkins & Binney 1 19901) would have 
blurred the kinematical borders of the two discs and finally 
merge them. Additionally, radial m igratio n due to transient spiral 
density waves (Sellwood & Binney'2002'), the interaction among 
multiple spirals (Minchev & Quillen 2006 |) or b ar and spirals 
(Minc hev & Famaev l20Tol Minchev et al. boilal Brunetti et al. 
120111). and the effect of mmor mergers (Quillen etal. l2009[ Bird 
et al. 201 ll) . introduces an even harder problem, since by chang- 
ing their angular momenta, stars arriving to the solar neighbor- 
hood have kinematics indistinguishable from those bom i n-situ. 
In case of satelli te accretion (e.g., Villalobos & Helmi l2008l 
Abadi et al. I2OOI the kinematical borders would also be blurred 
(because the consequent kinematical heating of the pre-existing 
disc) with, in addition, a chemical overlap of the accreted stellar 
populations to the Galactic population. High-resolution spectro- 
scopic surveys would be necessary to distinguish the chemical 
fingerprints of the extra-galactic population from the one born in 
situ. 

Both stellar heating and migration are time-dependent pro- 
cesses. Therefore, in the hypothesis that originally thin and thick 
disc were kinematically distinguishable, the natural question 
arises: How long does it take to delete the original kinematical 
borders between them? On the basis of our kinematic and chem- 
ical data, we infer that this time should be rather long (compara- 
ble to the Galaxy's lifetime), given that the RAVE sample con- 
tains stars with signs of kinematic diffusion together with the 
expected thin- and thick-disc populations. Models of our Galaxy 
togeth er with tools that create synthetic surveys (Sharma et al., 
I2OI ih will be employed in a next work in order to compare ob- 
servations with up-to-date models of the Milky Way. Qualitative 
and quantitative comparisons of our data with detailed chemo- 
dynamical models which follow the evolution of a Milky Way- 
like disc for the entire expected thin- and thick-disc lifetimes, 
are necessary to understand better the formation and evolution 
of the Galactic discs. 
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